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We investigate how the multiple bands of fermions on a crystal lattice evolve if a magnetic field is 
added which does not increase the number of bands. The kagome lattice is studied as generic example 
for a lattice with loops of three bonds. Finite Chern numbers occur as a nontrivial topological 
property in the presence of the magnetic field. The symmetries and periodicities as a function of 
the applied field are discussed. Strikingly, the dispersions of the edge states depend crucially on the 
precise shape of the boundary. This suggests that suitable design of the boundaries helps to tune 
physical properties which may even differ between upper and lower edges. Moreover, we suggest a 
promising gauge to realize this model in optical lattices. 


PACS numbers: 03.65.Vf,71.10.Fd,02.40.Pc,03.75.Lm 

I. INTRODUCTION 

Since the discovery of the quantum Hall effects mu, 
interest in topological aspects of condensed-matter sys¬ 
tems has risen and has stayed very high ever since. The 
link to topological invariants is given by the Berry phase 
[3] of the ground-state wave function in its dependence of 
magnetic fluxes on a torus [4 , 5 . This line of argument 
even shows that there are no corrections to Ohm’s law on 
the linear relationship between voltage and current mm- 

Only a few years ago, the discovery of topological in¬ 
sulators [8] gave another impetus to the field of topology 
in condensed-matter physics. Even in classical mechan¬ 
ical systems topological edges states can be excited [9j. 
Very recently, topologically nontrivial bands were real¬ 
ized in lattices of ultracold atoms mm and advocated 
in strongly frustrated spin systems with Dzyaloshinskii- 
Moriya interactions [T2| . 

For the present article, we are inspired by the realiza¬ 
tion of nontrivial gauge fields in optical lattices filled by 
ultracold atoms uni hb and by the interest in nontrivial 
lattices. Whereas Aidelsburger et al. realized bands with 
Chern numbers different from zero in Bravais lattices, we 
show that similar physics also occurs in crystal lattices, 
i.e., lattices with a basis. Moreover, our focus is on lat¬ 
tices with loops of odd numbers of bonds. The generic 
loop is a triangle which has the smallest odd number of 
bonds. In antiferromagnetic spin systems, it is the prime 
source of frustration. 

In our theoretical study, we aim at a proof-of-principle 
result. Though inspired by the impressive recent ad¬ 
vances in systems of ultracold atoms in optical lattices 
[lQl El El], we do not claim that this class of systems 
provides the most promising candidate for experimental 
realizations. The obstacle may be that the traps con¬ 
structed so far do not provide well-defined edges, which 
we examine below. But solid-state systems may fill the 
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gap. The recent observation of the quantum anomalous 
Hall effect in thin layers for ferromagnetic Chern insula¬ 
tors provides seminal progress mug. Presently, higher 
temperatures are reached at which the effect occurs m 
and theoretical calculations even suggest that in tailored 
systems experiments at room temperature are possible 
mm- There are concrete suggestions for tailored su¬ 
perlattices which should also make the design of partic¬ 
ular edges possible ii. An alternative to solid-state 
lattices may be the artificial lattices built from dots or 
antidots in tailored semiconductor structures [20]. 

In the present fundamental investigation, we study the 
kagome lattice shown in Fig. [l] This crystal lattice has a 
basis of three atoms, labeled A, B, and C in the sketch. 
One choice of primitive vectors is shown; they span a 
parallelogram which constitutes a unit cell. Note that 
the coordination number is only z = 4 which is rather 
low in comparison to z = 6 for the triangular lattice, the 
generic Bravais lattice with triangular loops. This lattice 
has appeared before in several studies. In particular, we 
show that the kagome lattice in certain magnetic fields 
corresponds to special cases of Haldane models [21], i.e., 
models without uniform magnetic field but complex hop¬ 
pings, on the kagome lattice [22H24] . Moreover, in the 
research field investigating fiat band models and inter¬ 
action effects on them the kagome lattice is also studied 
intensively; for a review see Ref. [25] . 

We study this model for the special magnetic fields 
where the number of bands still equals the number of 
bands without magnetic field. Note that this is an impor¬ 
tant difference from the investigation of the square lattice 
by Aidelsburger et al. where topological effects could only 
arise once the number of bands was increased by certain 
values of the magnetic field m- The bands and their 
dispersion are computed as well as the first Chern num¬ 
ber C n . This Chern number equals the Berry phase [3] 
occurring for a Bloch state which surrounds the Brillouin 
zone. The Brillouin zone is the relevant two-dimensional 
manifold, namely, a simple torus T 2 . This Chern number 
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FIG. 1. (Color online) Sketch of the kagome lattice, which 
is a crystal lattice with a basis of three sites (A, B, and C). 
Also shown is one choice of primitive vectors e* (in light blue) 
which span the unit cell. 

is defined by 

Cn := A /' F (n) df e z ( la ) 

Z7T1 Jji2 

— X 7 j ((^l^nk |^ 2 ^nk) (^ 2 ^nk \d\U n \^/^dk\dk2 • (lb) 

Z7T1 Jrp 2 

In Eq. ( [Ta| ) is the Berry curvature [3] of the princi¬ 
pal C7(l) fiber bundle defined by the nth Bloch eigenstate 
over the Brillouin zone. Equation ( |lb| ) is the explicit for¬ 
mula in terms of the Bloch states |n n k). The partial 
derivatives di refer to the derivation with respect to the 
momenta ki , i.e., di := d/dki. The above Chern num¬ 
ber must be an integer since it can be converted by the 
Stokes theorem to the integrated phase along a closed 
path divided by 2i r. 

We compute the above Chern number for the three 
kagome bands, thereby identifying the topologically non¬ 
trivial bands. 


II. MODEL AND MAGNETIC FIELDS 

For simplicity, we study a nearest-neighbor tight- 
binding model on the kagome lattice. Its Hamiltonian 
reads 

H = y^ tijc\Cj + ^2 V i°l C i ( 2 ) 

(ij) 1 

where the indices run over all sites and (ij) stands for 
nearest neighbors. The creation operator cj creates a 
fermion at site i and q annihilates it. In general, we con¬ 
sider complex hopping elements Uj = £exp(i$^) with 
some directed phases $ij resulting from the Peierls sub¬ 
stitution, i.e., 



where q is the charge of the hopping particle. The local 
potentials are given by Vi. They depend on the sublattice 
to which the site i belongs. 
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FIG. 2. (Color online) Dispersion of the three bands as func¬ 
tions of k\ and k<i as defined in Eq. 0 in the kagome lattice 
without any magnetic field and local potentials. The hopping 
is set to t — 1. Note the flat lowest band and the two Dirac 
cones where the two other bands touch with linear dispersion. 

As a reference, we first look at the case without any 
magnetic fields and without local potentials such that 
$ij = 0 and Vi = 0. It is known that in this case a flat 
band occurs because the kagome lattice is a line graph, 
namely, the line graph of the honeycomb lattice. Line 
graphs are generally known to have flat bands [26l|27]. 
The explicit calculation illustrates the flat band nicely in 
Fig- H It is based on the 3x3 matrix in momentum 
space 


( 0 cos&q cos& 2 \ 

COS&i 0 COS&3 J , (4) 

COS &2 COS &3 0 / 

where we use the wave-vector components 

ki := k • Aj, (5a) 

Ai := ei /2, (5b) 

A 2 := e 2 /2, (5c) 

A 3 := A 2 - A, (5d) 


with e 3 = (1,0) t and e 2 = (1, \/3)V 2 ! i- e -i the lattice 
constant from one A site to the nearest A site is set to 
unity. We stress that the full Brillouin zone is given by 
ki e [—7T/2, it/2} because the A* have only half the length 
of the primitive vectors. Note that k% = &2 — k\ holds 
according to the above definitions; only two of the ki are 
independent. 

The energy eigenvalues of the matrix 0 are 

Ti, 2 k = £(1 ± \A + 8 cos k\ cos &2 cos ks ), (6a) 

E 3k = -2 1. (6b) 

These bands are degenerate at k = (0, 0)^ where the flat 
band touches the lower dispersive band and at k± = 
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±(27r/3)(l, — \/3)t where the two dispersive bands dis¬ 
play Dirac cones. Due to the degeneracy, the Chern num¬ 
bers C n of the bands are not properly defined. But even 
if one breaks the degeneracy, for instance by infinitesimal 
local potentials, the C n turn out to be trivial, i.e., zero, 
because the model without finite phases is time-reversal 
invariant which implies C n = 0 . 

Thus we turn to the model in a uniform magnetic field 
perpendicular to the kagome plane. This induces finite 
phases according to Eq. 0 - We do not aim at discussing 
the complete intricate interplay of discrete lattice sym¬ 
metry and noncommuting magnetic translations [28b3Q] 
which give rise to multiple Hofstadter bands [31]. In¬ 
stead, we aim at the situation where the magnetic trans¬ 
lations Ti along the primitive vectors commute. They 
obey the relation T{I 2 = exp(i 0 uc)^ 2 ^i where 0uc is 
the magnetic flux through the unit cell measured in units 
of h/e [28ff3Qj . Thus, in order not to increase the number 
of bands, we avoid decreasing the translational symme¬ 
try. Thus we focus on the commuting case X 1 X 2 = X 2 X 1 
implying 


<t> uc = 27T n, (7) 

where n G Z is an arbitrary integer. Inspecting the unit 
cell shown by the lightly shaded area in Fig. [3] element ary 
geometry shows that 0 uc = 80a if 0 a is the flux through 
a triangle. Thus Eq. 0 translates to 

= jn. (8) 

The physics modulo gauge transformations of the phases 
is only influenced by the fluxes through the loops mod¬ 
ulo flux quanta, not by the individual phases. Thus we 
conclude that the physics, such as Chern numbers of the 
bands, depends on 0 a with periodicity of 27r so that we 
will have to study only the values of n modulo 8 . 


this point explicit we consider a gauge for arbitrary flux 
in AppendixjA] and show that it can be regauged to the 
one in Fig. [3] in Appendix [B] In this particular gauge, 
the relation between these phases and the flux 0 a is the 
following 


3di = 0 a — 27rrai, 
30 2 = 0A — 27T 77l 2 , 
-3($i + $ 2 ) 60a - 27T ra 3 , 

0 = 80a — 27rn, 


(9a) 

(9b) 

(9c) 

(9d) 


where rrii and n £ Z may occur since the fluxes through 
the loops are only fixed up to multiples of 27r. The last 
equation [Eq.(|9d|)] is the sum of the three equations be¬ 
fore with n = rrii + ra 2 + m 3 which confirms Eq. 0 . 
The left-hand side of Eq. (9d) stands for the vanishing 
sum of all phases around a unit cell; see Fig. [3] This 
does not mean that there is no net uniform magnetic 
field through the lattice, but it reflects the fact that the 
magnetic translations commute for the special magnetic 
fields we are considering because the flux through the 
unit cell is a multiple of the flux quantum. Thus, there 
are gauges without net total phase around the unit cells 
and we are employing such a gauge here for simplicity. 
But in the smaller loops, i.e., the triangles, there is a net 
flux given by 0 a- 

For finite phases and local potentials we obtain the 
following matrix problem: 


2 te“ i % 2 te^c 2 \ 
c 3 

c 2 2 te lu c 3 


V A 

H* = | 2te w c\ V B 2te~ w c 3 I , 


\2te ^c 2 2te 1 ' ,> cp, Vc ) 


( 10 ) 


where we use the shorthand 


Ci := cos (ki — (—1)*A0) 


(ii) 


B $2 A 


V ■ V ■ \A X/: W 

aa\ . aaa\: :. a 

A 0! B 0 X 0 X #1 e i 01 


FIG. 3. (Color online) Gauge chosen for the specific calcu¬ 
lations. Two phases, 0i and 02, are introduced which are 
directed in the sense of the arrows. 

We choose a specific gauge where the phases 0^- take 
values as depicted in Fig. [3] The arrows on the bonds 
stand for the direction of the phases 0i and 0 2 . In this 
gauge, the total phase around the unit cell takes the value 
zero, which may be a surprise because it seems to reflect 
zero flux and thus zero magnetic field. But this is not 
the case because the sum of all phases along closed loops 
is meaningful only modulo 27r. Hence zero total phase 
complies with a finite flux fulfilling Eq. ([ 7 ]). To make 


for brevity. We have introduced the average and the 
difference of the phases 0 := (0i +02)/2 and A0 := (0i — 
0 2 )/2. Obviously, the matrix Hk is 27r periodic in 0 and 
A0. Interestingly, the band Hamiltonian (10) is identical 
to the ones considered in Haldane models before j22lf24j if 
A0 = 0 and the local potentials are switched off. But we 
recall that Eq. (10) only holds for uniform magnetic fields 
which meet the condition (|8|. This observation leads us 
to the experimentally useful result that certain Haldane 
models can be realized for specific fluxes complying with 
Eq. |8| by uniform magnetic fields. 

But the eigenvalues are even periodic in 0 with period 
27 t/ 3. This can be seen by the gauge transformation 


(10 0 \ 

U = 0 e i27r / 3 0 , (12) 

\0 0 e i47r / 3 / 


which transforms the Hamiltonian matrix according to 
740? + 2 tt/3) = UH\ (13) 
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The value of A$ remains unchanged. Clearly, the eigen¬ 
values at each momentum do not change under the trans¬ 
formation and hence they are 27r/3-periodic in$. Since 
the transformation U is independent of momentum it 
does not change the Chern numbers either so that the 
Chern numbers are also 27r/3-periodic. Note that a 
change of & by 2ir/3 at constant corresponds to the 
simultaneous change 


m i 


■ mi— 1,7712—>> m 2 —1, and m 3 m 3 + 2. (14) 


Alternatively, the increase of $ by 2tt/3 can be inter¬ 
preted as an increase of 0 a by 2i r and the increment 
m 3 —>• m 2 + 8 . Thus, we see that the physics is indeed 
27r-periodic in 0 a- 

In addition, one easily sees that $ $ + 7 r simply 

leads to a global minus sign in the Hamiltonian matrix 
( 10 ), thus inverting the sequence of bands and thereby 


swapping the Chern numbers of the first and the third 
bands. The change & $ + 7r directly corresponds to the 

change 0 a —» 0a + 37 t; see Eq. ©• But in view of the 
27r-periodicity in 0 a we find the corresponding bands up 
to a gauge transformation already for 0 a 0a + tt. For 
examples, we refer the reader to the Sec. |III| below. 

Moreover, reversing time amounts up to changing the 
sign of the phases tk —>• —tk and of the flux 0 a — 0 a; 
i.e., the Hamiltonian matrix is complex conjugated. This 
transformation of Hk leaves the eigenvalues and, thus 
the dispersion unchanged because the energy eigenvalues 
are real numbers. But it changes the sign of the Chern 
numbers. 

Combined with the shift of the flux by 7 r we obtain the 
result that the change 0 a ^ tt — 0 a inverts the energy 
bands, i.e., takes the dispersion energy to their negative 
values, but leaves the Chern numbers unchanged. We 
will illustrate these properties in Sec. m below. 

The Chern numbers of the bands are also 27r/3-periodic 
in A'd because the shift of the momenta 


h\ —y k\ T 27r/3, 
k 2 k 2 — 27r/3, 


(15a) 

(15b) 


implies k% —» ks — 47r/3 = k% + 27r/3, where the last iden¬ 
tity holds modulo 27r. Thus, this shift of the momenta 
corresponds to the change of gauge Ad —»• A'd + 27r/3. So 
the eigenvalues at given momenta change, but they are 
only shifted in reciprocal space. Since the Chern number 
is a property of the band in the Brillouin zone it does not 
change under Eqs. (15). This is obvious if one uses the 


shifted Brillouin zone, shifted by the same amount as in 
Eqs. (15). Note that a change of Ad by 2tt/3 at constant 
d corresponds to the simultaneous change 


mi 


mi + 1 , rn 2 —>> m 2 — 1 , and m 3 —» m 3 . (16) 


We stress that the combination of the two gauge trans¬ 
formations, Eqs. (14) and (16), illustrates that the essen¬ 


tial physics is determined by n = mi + m 2 + m 3 while the 
individual values of the m^ do no matter much. However, 
m 3 cannot be regauged to change by unity only. 


III. RESULTS 

In this section, we show explicit results for the disper¬ 
sion and the Chern numbers of the three bands of the 
kagome lattice in a magnetic field obeying relation (| 8 |. 


A. Results without local potentials 


For vanishing local potential the eigenvalues of the ma¬ 
trix (10) can be determined analytically. Solving the 
characteristic polynomial of order 3 we obtain 


E bk = 4i y/Q cos ([0k + 2 tt5]/3) , b e {0,1, 2 }, (17a) 

0k := arg (Fcos(3$) + i03 3 — [Pcos(3i9)] 2 ^) , (17b) 

P := cic 2 c 3 , (17c) 

Q := (c? +4 + c§)/3. (17d) 



FIG. 4. (Color online) Dispersion of the three bands as func¬ 
tions of ki and k 2 as defined in Eq. (pi in the kagome lattice 
at 0 a = 7r, i.e., n — 4 in Eq. ([8]), and tk = 0a/3 and no local 
potentials. The hopping is set to t = 1. 


For simplicity, we choose mi = m 2 = 0, thus Ad 
0. and vary m 3 = n to achieve various fluxes 0a- As 
pointed out above, one only has to consider eight values 
of 0 a because of the 27r-periodicity and the quantization 
condition (| 8 |. Additionally, 0 a —^ 0a+tt only inverts the 
energies and swaps Chern numbers so that it is sufficient 
to consider n = 0,1, 2, 3. For instance, n = 4 corresponds 
to adding the phase 7 r and, indeed, Fig. [4] displays the 
negative bands of Fig.|2j The degenerate, touching bands 
do not allow for an unambiguous definition of the Chern 
number. But the fact that the Hamiltonian matrix is 
equivalent to a real matrix up to a gauge transformation 
(14) teaches us that the bands are topologically trivial 
C n = 0. 

The three cases n = 1,2,3 are much more interesting; 
the corresponding bands are depicted in Fig. [5] in ascend¬ 
ing order. All three bands are separated from one another 
so that the Chern numbers are well defined. Note the 
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FIG. 5. (Color online) Dispersion of the three bands as func¬ 
tions of ki and k <2 as dehned in Eq. (pi in the kagome lattice at 
0 a — (7r/4)n for n = 1 (top), n = 2 (middle), and n — 3 (bot¬ 
tom) and — 0a/ 3 and no local potentials. The hopping is 
set to t = 1. 


gradual evolution upon increasing flux, illustrated nicely 
from Fig. [2] over the panels of Fig. [5] to Fig. [4j 

The fact that 0 a tt — 0a inverts the signs of the 
energies is nicely illustrated by comparing Fig. [2] and Fig. 
[4j The same relation is seen between the upper and the 
lower panels of Fig. [5| The middle panel of Fig. [5] remains 
unchanged under inversion of the signs, which implies 
that the upper and lower bands differ only by their sign 
and that the middle band is identical to zero. 

Next, we turn to the Chern numbers of the bands. We 
did not find a way to evaluate the expressions in Eqs. ([l]) 
analytically. Moreover, the two-dimensional integrals are 
difficult to implement numerically to high precision. But 
due to the robustness of the discrete Chern numbers the 


computation of accumulated phases on a finite mesh is 
very accurate, even if the mesh is not particularly dense 
{32j. Thus we use this approach to compute the Chern 
numbers shown in Fig. [6j 



FIG. 6. (Color online) Chern number C n as defined in Eqs. 
0 computed for various fluxes 0 a for = $2 = 0a/ 3. The 
bands are numbered in order of ascending energy. 

Several interesting observations can be made. First, 
the fluxes, which are multiples of 7 r, indeed yield trivial 
bands. Second, incrementing the flux by 1 r yields a swap 
of the Chern numbers due to the inverted global sign of 
the energies. The change 0 a —^ tt — 0a leaves the Chern 
numbers unchanged. Third, swapping the sign of the 
fluxes swaps also the sign of the Chern numbers. Fourth, 
the gradual evolution of the bands from 0 a = 7r/4 to 
0a = 37t/4 with energetically well-separated bands is 
reflected in constant Chern numbers. 

Further values of 0 a do not need to be analyzed be¬ 
cause of the 27r-periodicity of the Chern numbers as a 
function of 0 a- Thus the pattern in Fig. [6] will occur 
repeatedly. 


B. Robustness against local potentials 

In order to illustrate that the Chern numbers are ro¬ 
bust against small changes of the Hamiltonian we turn 
to the case with local potentials. Clearly, equal changes 
of all potentials Va,Vb , and Vc do not affect the Chern 
number at all because they simply add an identity ma¬ 
trix to Hk which shifts the eigenenergy but has no impact 
on the eigenvectors and hence no impact on the Chern 
numbers; see Eqs. 0 . 

For illustration, we study two nontrivial patterns of 
the local potentials, namely Va = S = —Vb and Vc = 0 
in Fig. [7] and Va = 5 = Vc and Vb = 0 in Fig. [8] These 
patterns are applied to the case 0a = tt/2 of which the 
dispersion without local potentials is shown in the middle 
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panel of Fig. [5] Clearly, the nontrivial Chern numbers of 
the lower and upper bands are robust against the local 
potentials. It is necessary to apply potentials of the order 
of S ~ t to destroy the topological phases because these 
phases are protected by gaps; see the middle panel of Fig. 
[5] This is numerically shown in Figs. [7] and [5] and can be 
analytically verified. For the pattern Va = S = —Vb and 
Vc — 0, the gaps close at |<S| = and for the pattern 
Va = S = Vc and Vb = 0, they close at \S\ = 2. 

These exemplary results illustrate the robustness of 
topological phases against perturbations. For the topo¬ 
logical character to be lost the bands must touch and be¬ 
come degenerate. A deformation alone is not sufficient. 



FIG. 7. (Color online) Chern number C n as defined in Eqs. 
0 computed for flux 0 a = 7r/2 for #1 = $2 = 7r/6 at finite 
local potential according to Va = S = — Vb and Vc — 0. The 
hopping t is set to 1 as before. 



FIG. 8. (Color online) Chern number C n as defined in Eqs. 
0 computed for flux 0 a = 7r/2 for 'di ='$2 = 7r/6at finite 
local potential according to Va = S = Vc and Vb — 0. The 
hopping t is set to 1 as before. 


IV. EDGE STATES 

For measurable quantities the bulk properties are not 
of prime interest. The main difference between a usual, 
trivial band insulator and a topological one occurs at 
the boundaries where states with vanishing gap have to 
appear since otherwise the integer Chern numbers cannot 
change from their finite value to zero. These edge states 
are crucial for the interesting properties of topologically 
ordered systems. 



FIG. 9. (Color online) Three types of edges: with triangular 
teeth at both edges (top), for smooth boundaries without tri¬ 
angular teeth at both edges (middle), and for mixed bound¬ 
aries which are smooth at the upper edge and saw-toothed 
at the lower edge (bottom). The bracket labeled U N units” 
stands for N times repeated units between the boundaries. 

To illustrate the existence of such edge states also in 
the kagome lattice in a magnetic field we diagonalize 
strips of it. A strip is infinitely extended in the x direc¬ 
tion such that k x is a conserved quantum numbers. But 
in the y direction the strip is of finite extension. Three 
different situations of the boundaries are displayed in Fig. 
[9] The height of the strips is fixed by the number N of re¬ 
peated units of six sites in the perpendicular y direction. 
The dispersions of the in-gap edges states converge very 
quickly for N —» 00. Our results are based on computa¬ 
tions with N = 80 which amounts to the diagonalization 
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FIG. 10. (Color online) Dispersions of the edge states at 
0 a = 7r/2 if the edges are shaped as shown in Fig. [ 9 ] the top, 
middle, and bottom panels here shows the data relevant for 
the corresponding panel in Fig. [9] The right-moving states 
from the lower edge are shown in orange and the left-moving 
ones from the upper edge in dark green. The shaded area 
stands for the continua of the bulk states if k y is not fixed. 
The dashed line at zero stands for the nondispersive bulk 
states of the trivial middle band; cf. middle panel in Fig. [5] 


of matrices of dimension 480. 

The resulting dispersions are shown in Fig. 10 Strik¬ 
ingly, there occur significant differences for the disper¬ 
sions of the differently shaped boundaries. Even the 
group velocities differ by about a factor of 2. 

These differences make the assignment of the edge 
states to the upper and lower edges particularly simple 
by analyzing mixed boundaries; see the lower panels in 
Figs.[9]and[ToJ Clearly, we partly find the dispersion from 


the upper panel for the right-moving states and from the 
lower panel for the left-moving states. Thus the right- 
moving states live at the lower edge and the left-moving 
states at the upper edge. Of course, this will be swapped 
if the sign of the magnetic field is inverted. 

We conclude that the shapes of the boundaries open an 
interesting field for tuning the properties, for instance the 
transport properties, of topologically nontrivial bands. 


V. CONCLUSIONS 

We have shown that a crystal lattice with basis equally 
allows us to induce topologically nontrivial bands upon 
application of a uniform magnetic field. Interestingly, the 
number of bands does not need to be increased from the 
values at zero magnetic field, which is in contrast to the 
situation in Bravais lattices HU Ell- 


Y'Y Y'Y. Y 

A. A. A. A:a 


FIG. 11. (Color online) Phases after the gauge transforma¬ 
tion c\ c\ exp(=pi$) for A and B sites, respectively. The 
experimental realization should be easier in this gauge. 

In particular, we studied nearest-neighbor hopping on 
the kagome lattice with a three-fold basis. This lattice 
is made of triangular loops, i.e., loops with the smallest 
number of odd bonds which generically induces frustra¬ 
tion effects. Without magnetic field the bands are not 
separated, but the special magnetic fields which preserve 
the commutation of the translations induce a proper 
splitting and nonzero Chern numbers for two of the three 
bands. Interestingly, we found that the kagome lattice in 
a uniform magnetic field of particular strengths realizes 
Haldane models on this lattice for particular quantized 
fluxes. 

It would be interesting to realize crystal lattices with 
nontrivial gauges experimentally mm- In order not 
to be forced to induce finite phases P := #1 = $2 for 
all hopping processes, we point out that a simple gauge 
transformation eliminates all phases on bonds to or from 
sites of the sublattice C. It consists of the additional fac¬ 
tor ct — »• ct exp(— ifi) for A sites and the additional factor 
ct —>• ct exp(iP) for B sites. Of course, the fluxes through 
the triangles are left unchanged because the phases on 
bonds between the A and B sites are tripled, $ —>> 3P, as 
shown in Fig. El 

Finally, we investigated the edge states for three ways 
to cut the kagome lattice along the same direction. 
Again, this possibility is a particular feature for lattices 
with a basis. Surprisingly, we found a significant quali¬ 
tative dependence of the dispersion of the edge states on 
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the nature of the boundaries. The shape and the position 
of the dispersions differ as well as the group velocities. 
Still there are left- and right-moving modes which can be 
clearly assigned to the upper or lower edge. The situation 
is particularly clear if the upper and lower edges differ. 
Then we obtained two differing edge states, one reflecting 
the behavior at the upper edge and one the behavior at 
the lower edge. 

Experimental verification of this interesting scenario 
is called for. One may envisage fascinating applications 
because the transport properties will differ depending on 
the direction of motion of a fermion from left to right or 
vice versa. It is conceivable to build diode like devices if 
the properties of the edges are tuned to maximize their 
differing characteristics, for instance the group velocity. 

At present, the most promising systems for experimen¬ 
tal realization are not yet clear [25] . But recent advances 
in realizing artificial gauges in systems of optical lattices 
mamma make it likely that soon the theoretical results 
presented here can be put to experimental tests. The re¬ 
alization of precisely defined edges may pose a problem, 
but any reproducible difference between two edges of a 
strip of lattice will make a verification possible. 

Alternatively, various solid-state systems may consti¬ 
tute alternative routes towards experimental realizations, 
for instance thin slabs of ferromagnetic Chern insulators 
mm, honeycomb lattices built from atoms with higher 
atomic number na. artificial lattices on surfaces mm, 
or special lattices built from semiconductor nanostruc¬ 
tures EDI . In any case, we conclude that designed crystal 
lattices constitute a promising field of research. 
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Appendix A: Arbitrary magnetic flux 


In Fig. 12 a general gauge is shown which is a lattice 
version of the Landau gauge. The phases depend on the 
x coordinate of their bond, i.e., on the x coordinate of 
the midpoint of their bond. The dependence is linear, 


ip{x) := 40a^, 


(Al) 


where we imply that the lattice constant is set to unity. 

Then it is easy to see that indeed the flux through a 
triangle is given by 


<t >A = <P{x + -) - <p(x) 


(A2) 


so that our notation is consistent. Similarly, we find for 



FIG. 12. (Color online) Gauge for arbitrary magnetic flux, 
i.e., for arbitrary homogeneous magnetic field. Recall that the 
lattice constant, i.e., double the distance between neighboring 
sites, is set to unity. 


the flux through the hexagon 


0hexagon — 2(p(x + 1) 2(p(x T ^) (A3a) 


= 60 a 


(A3b) 


as it has to be for a homogeneous magnetic field. In 
total, the flux through a unit cell is given by 80 a- Note 
that these results do not change if we consider shifted 
triangles or hexagons in other parts of the lat tice . Shifts 
in the y direction do not change the phases ( |A1| ) at all. 
Shifts in the x direction do not change the fluxes because 
they depend on phase differe nces being independent of 
x due to the linearity of Eq. (Al). So the gauge shown 


in Fig. [12] correctly describes the effect of an arbitrary 
homogeneous magnetic field perpendicular to the plane 
of the lattice. 


Appendix B: Translational invariant gauge for 
specific magnetic fluxes 


Here we show that for 0 a = n7r/4 the gauge in Fig. 12 


can be regauged to yield the gauge in Fig. [13} For clarity 


we focus on the case $1 = $2 = 0a/3, i.e., mi = m 2 = 0 
in Eq. (|9|. 

We leave the fermion at site a in Fig. [13] unchanged, 
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FIG. 13. (Color online) Gauge for specific magnetic fluxes, 
i.e., for specific homogeneous magnetic fields. The necessary 
regauging on the various sites a to h and a' to e are given in 
the main text. 


but we transform the others according to 


c ^ — y 

(Bla) 

c l c t ex p( il? i)> 

(Bib) 

c\ —y <4exp(2i0i - iip(x + 1/4)), 

(Blc) 

c\ —y 4exp(3i0i - i2ip{x + 1/4)), 

(Bid) 

4 —y 4 exp(i0i — icp(x) — i(p(x + 1/4)). 

(Ble) 


The annihilation operators are transformed by the com¬ 
plex conjugate phase factors. Obviously, the phase fac¬ 
tors are chosen such that the phases on the bonds <26, 6c, 
ce, and cd match those in Fig. [13] We have to check the 


resulting phase factor on ca and compute 

d ca = -2#i + <p(x + 1/4) - <p(x) (B2a) 

= -(2/3)0 a + 0a (B2b) 

= 0 a/3 (B2c) 

= ^i, (B2d) 

which is what we wanted to have. Analogously, we con¬ 

firm for the phase on bond ed 

$ed = &i ~ <p(x) - 3#i + ip(x + 1/4) (B3a) 

= -(2/3)0 a + 0a (B3b) 

= 0a/3 (B3c) 

= tf 2 . (B3d) 


If we now pass from the sites a to e to the site s a' to e! 
we can use almost the same transformations (B1) because 
the phases on the bonds are changed by multiples of i r 
due to 


(j)(x + 1) = (j>{x) + 40a (B4a) 

= (j>{x) + nir (B4b) 


for any x. This is so because 0 a takes the particular val¬ 
ues mr/4 representing the crucial step in our argument. 


If n is even we may use exactly the same transfor¬ 
mations. If n is odd we regauge the center site c' by the 
factor — 1 and we are back to the phases between the sites 
a through e and use the regauge transformation (Bl). In 
addition, it is obvious that the phases on the bonds bo! 


and ed! correspond to the ones in Fig. 13 


For completeness, we state the regauge transforma¬ 
tions on sites /, g, and ft, 


cj. -T cj- exp(2i$i — 3i(p(x) — 4z0a), (B5a) 

c l c l ex P(^i — 4 iif{x) — 6i0 a), (B5b) 

4 -T c t exp(3i0i — 4 icp(x) — 7z0a), (B5c) 


because they are not directly deduced from the ones on 
a to e. Straightforward calculations verify that these 
transforms yield the phase pattern shown in Fig. |T3j The 
sites resulting from a shift by one lattice constant to the 
right are again transformed either exactly the same way 
as the sites e to h and d! if n is even. For odd n, the 
same transformations are used except for an additional 
sign change on the shifted center site /'. 


In the next upper row the phase pattern for the row of 
sites a to e is repeated. Thus one can reuse the transfor¬ 
mations (Bl) except that all the phas es mu st be taken 
relative to the phase of the h site [Eq. (B5c)]. 


In an analogous fashion, one can implement regauge 
transformations yielding ^ $ 2 , but they are slightly 
more complex because they comprise a certain phase shift 
upon x x + 1. 


In conclusion, we can regauge the case of a general 
homogeneous magnetic field to the phases shown in Figs. 

|3]or[Il 
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